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Neutrinos propagating tlirougli matter may participate in forward coherent neutral-current-like 
scattering arising from non-standard interactions as well as from the Mikheyev-Smirnov-Wolfenstein 
matter potential Ve. We show that at fixed long baselines through matter of constant density, the 
non-standard interaction potential e^r can contribute an additional term to the oscillation phase 
whose sign differs for versus propagation in matter. Its presence can cause different apparent 
Am'^ to be erroneously inferred on the basis of oscillations in vacuum, with values lying above (for 
Vfj) or below (for v^i) the actual /^m%2 for the case where e^r is predominantly real- valued and of sign 
opposite to Am|2- An NSI scenario invoking only 5R(e^T) is shown to be capable of accounting for a 
disparity recently reported between oscillation survival for 77^ and i^^ fluxes measured at 735 km by 
the MINOS experiment. Implications for mantle traversal by atmospheric neutrinos are examined. 
The NSI matter potential with non-maximal mixing could evade conventional atmospheric neutrino 
analyses which do not distinguish from on an event-by-event basis. 
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I. INTRODUCTION 

Accelerator-based neutrino long baseline experiments 
have recently entered an era wherein oscillation mea- 
surements using antineutrino exposures are being under- 
taken. Since previous experimental investigations have 
been have based in the main on neutrino beam exposures, 
the on-going and next-round explorations using antineu- 
trino reactions will provide a crucial complementary per- 
spective on oscillation phenomena. The question natu- 
rally arises as to whether standard phenomenology, de- 
scribing mixing among either two or three of active flavor 
neutrinos and without or with inclusion of the Mikheyev- 
Smirnov-Wolfenstein (MSW) matter effect provides 
an adequate framework - or whether a more expansive 
formalism is required. In this work we investigate a mat- 
ter effect scenario which is based upon a specific neu- 
trino non-standard interaction (NSI) which we identify 
as highly promising for generating distinctive observable 
effects in the neutrino "atmospheric" /i-r flavor mixing 
sector. 

Impetus for this paper arises in part from observa- 
tions of muon-neutrino flavor disappearance for 77^ — > 77^ 
and for — >■ obtained by the MINOS experiment. 
As recently reported 0, MINOS infers the value for 
Am|2 governing oscillations at 735 km to be some- 
what higher than the value gleaned from disappear- 
ance. The new observations are in line with the MI- 
NOS low statistics result of last year, in which 
survival was examined using antineutrino (background) 
reactions isolated in exposures with neutrino beams Q. 
While the earlier observation was of poor statistical sig- 
nificance, it attracted attention as a possible harbinger 
of new physics, e.g. of CPT violation in the neutrino sec- 
tor [3] , or of new interactions with matter-driven mixing 
to a sterile neutrino [1], or of an NSI matter effect 0. 

We begin with a brief summary in Sec. II of salient 
points from the literature on NSI matter effects as they 



pertain to neutrino oscillations. In Sec. Ill we then spec- 
ify the neutrino oscillation phenomenology which follows 
from existence of a single NSI matter potential e^rVe, 
proposed to exist in addition to the MSW matter effect. 
The latter, well-known matter potential is Ve = y/2Gpn^, 
where Gi? is the Fermi constant, and Ue is the electron 
number density averaged over the neutrino path through 
the Earth. This sets the stage for discussion of certain 
remarkable implications for experiments in Sec. IV. Our 
formalism applies most directly to neutrino traversal of 
the Earth's crust, and we use recent results from the MI- 
NOS long baseline experiment to illustrate such an ap- 
plication. We proceed to show that a real-valued e^^ in 
conjunction with a single Am^2 ^^'^ with a single mix- 
ing angle 623 provides a satisfactory description for the 
MINOS antineutrino and neutrino beam data currently 
available. In Sec. V we relate our formalism to neutrino 
traversal of the Earth's mantle, and identify effects of 
relevance to measurements using atmospheric neutrinos. 
Section VI concludes the paper with mention of near- 
future experimental developments which can shed light 
on the existence of e^,- , and on other new physics mecha- 
nisms which also have capability to differentiate from 
//^ oscillations. 

II. NEUTRINO NSI MATTER EFFECT 

Under the hypothesis that neutrinos participate in 
heretofore unobserved non-standard interactions, the 
neutrino flavor Hamiltonian will carry new potential 
terms analogous to the MSW potential. The latter poten- 
tial accounts for coherent forward scattering of electron- 
flavor neutrinos from the electrons within a field of mat- 
ter. Possibilities for the underlying neutrino scatter- 
ing processes include both flavor-changing and flavor- 
conserving NSI. There are six possible NSI amplitudes 
which can arise in neutrino propagation through matter, 
conventionally designated as Cee, ^e^i, Cer, f/^r, and 



2 



Stt- CP-violating effects are possible in the event that 
the amplitudes carry phases. The effects which may arise 
from various combinations of eap have received extensive 
treatment in the literature; recent summaries with rele- 
vant references can be found in Refs. 0j0|- 

For oscillation survival probabilities of the (f^, i^r) sec- 
tor of interest here, direct contributions may arise from 

the three potentials e^^K, e^^Ve, and e^^Ve- In a i^^ on an event-by-event basis is available, there is a risk 



ing limits, while possibly disfavored, are not as yet ruled 
out. 

As we elaborate in Sec. V, the phenomenology allows 
structure which only becomes apparent when sam- 
ples are treated separately from i/^ samples 16]. More 
generally, our study suggests that, in fitting an NSI e^,- 
scenario to data for which no distinction between and 



recent evaluation of model-independent bounds on pro- 
duction and detection of neutrino NSI it was found that 

fj,r\ < 0.33 and Ett is even less 



< 0.064, whereas |e 
constrained 



We note that limits on |e^^| and |e^^| 



of "averaging" over NSI structure in such way as to re- 
duce to a mimicry of conventional vacuum oscillations. 

III. TWO-FLAVOR EVOLUTION WITH 
MATTER EFFECT 



To obtain accessible expressions of sufficient accuracy, 
we neglect solar-scale mixing {Ani^i = 0) and subdom- 
inant z/^ — ^ Ve oscillations (^13 = 0), and work in a 
two-state mixing framework. Unless otherwise noted, 
the normal hierarchy for neutrino mass eigenstates is 
assumed. Basis states for neutrino flavor (a = /x, r) 
are then , and the evolution of states in time 
(h ^ c ~ 1, so t = L for ultrarelativistic neutrinos) is 
governed by the effective wave equation 



based on solar and KaniLAND neutrino data are re- 
ported from an analysis which restricts to neutrino NSI 
with electrons (but not with u or d quarks) [lol |. 

In this work we neglect the flavor-conserving NSI and 
focus on e^T-K as the single NSI capable of producing 
significant differences in versus Vij, oscillations in long- 
baseline experiments. The extraordinary role that e^r 
can play in distinguishing from oscillations has re- 
ceived little direct discussion in the literature to date. 
This circumstance reflects in part a dearth of sufficient 
experimental information to motivate NSI expositions to 
get down to specifics. However a recent work has dis- 
cussed a somewhat analogous role for the Cer NSI (of 
magnitude ~10%) with explaining discordant 0i2 val- 
ues inferred from solar (electron) neutrinos versus Kam- 
LAND (electron) antineutrinos 11]. 

A two-flavor neutrino mixing framework is an adequate 
venue for examination of and propagation in the 
constant-density terrestrial crust, and we invoke such a 
framework in Sec. III. It has been argued, on the ba- 
sis of fitting using the two-flavor mixing framework to 
SuperKamiokande and MACRO atmospheric (z^^ + T^fi) 
data, that the trend of muon-flavor survival at high ener- 
gies and for global baselines constrains \ef_,r\ to be a few The full Hamiltonian is then given by 



(1) 



The Hamiltonian in flavor basis describing neutrino 
propagation in vacuum is obtained from the vacuum 
Hamiltonian in mass basis H^^^^ = diag(0, Aml2/2E^) 
through a rotation via the standard unitary mixing ma- 
trix i?i (023 ) . We augment the flavor basis vacuum Hamil- 
tonian with the NSI potential term 
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e^.K 
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percent or less [12|, |13[. For neutrino and antineutrino 
propagation over globe-spanning baselines however, there 
are complications. Conventional three-flavor oscillations 
requires treatment of the different matter densities of the 
crust, mantle, and core, with the MSW resonances of 
the mantle and core coming into play. With inclusion 
of NSI matter effects into the "mix", e.g. Cer as well 
as e^T, significant new degrees of freedom become avail- 
able which can couple to ones neglected in a two-fiavor 
mixing framework. It has been demonstrated for eer, 
Ctt, and Cee, that NSI bounds derived from atmospheric 
neutrinos via two-flavor analysis become relaxed when Vf, 
mixing is included [l3, llEl ■ The possibility also remains 
that a similar outcome would ensue in a full three-flavor 
analysis extending to e^r- Consequently, the verity of 
stringent e^r bounds obtained using the two-flavor mix- 
ing framework for atmospheric neutrinos has been called 
into question . A three-flavor framework treatment 

for atmospheric neutrinos with NSI is generally regarded 
to be preferable [l^. In the absence of a comprehen- 
sive three-flavor framework analysis at this time, we take 
the view that le^rl values exceeding the two-flavor mix- 



(2) 



(3) 



Only the real part of e^r distinguishes between and 
F^ in the derivations to follow. Hereafter we neglect the 
CP-violating imaginary part of e^r and assume e^r = 
3fi(e^r) in order to focus on the physics implied by the 
real part of e^r ■ 

After algebraic manipulation leading to removal of a 
term proportional to I, which merely contributes an over- 
all phase to the oscillation amplitudes, the Hamiltonian 
can be expressed as 



where a is the Pauli vector, and 



(4) 



'32 



-t-e^^V;,0,-cos(26'23)- 



'32 



(5) 

The corresponding evolution operator U{t — L,0) is 
U{t, 0) = e"*^*"'* = e-'^-^ ^ = e"*"-'^ . (6) 
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The evolution operator U amounts to a rotation of 20 
about the direction n = iV/|A^| in the two-flavor spinor 
space, where (j) = \N\L \s given by 



'32 



T2sin(2023) 



Am 



1/2 



(PoT) and 7.2 x 10^" PoT using Reverse and Forward 
Horn Current running (hereafter referred to as the RHC 
and FHC exposures). It is reported that the mass- 
squared difference of neutrino mass states, Arn^ , deduced 
from the observed survival rate at the MINOS far de- 
tector at 735 km has best fit value of 3.36to;4oX 10~^ eV^ 
This value is higher than the updated MINOS best fit 



value for disappearance. Am =2.35 



+0.11 
-0.08 



X 10--^ eV^ 



+ {e,rV,Y\ L. (7) 

Both j/^ and data are compatible with maximal mix 



ing, however the best fit for z/^ allows more deviation: 



In Eq. © the middle term admits both signs, depend- sin" 26 = 0.86 ± 0.11 for i/^ compared to sin 29 > 0.91 



ing on whether it describes z^^ or i^^. This is because 
the matter potential becomes negative for 17^ propa- 
gation. (In Eq. ©, and in Eqs. ©, (HU]) to follow, the 
upper (lower) sign refers to (t'p) oscillations, whereas 
Cfj^r carries its own fixed sign.) 

From the time evolution operator, the survival proba- 
bility of (Fp) at i = L is 



^ (77) 



1 -J"sin^ 



where 



cos^ 261 



23 



l + 2sin(2M%^i?, 



(8) 



(9) 



It can be seen for the case e^r — that reduces to 
sin^ 2023, and (p reduces to the vacuum phase, hence the 
vacuum survival probability is recovered. For the case 
of maximal mixing, 023 — 45° and cos 2023 = 0, the 
factor becomes 1 and the survival probabilities simplify 
to 



sm 



Am32 _ 



4E^ 



eur\Ve\ 



L 



We observe that relative to the standard survival proba- 
bility, the vacuum oscillation phase is augmented by an 
additional term which is independent of E^, ; it is a mat- 
ter effect induced by the flavor-changing NSI. It is impor- 
tant that the NSI matter effect phase term appears with 
a different sign in the expressions for and 77^ survival 
probabilities. 

Our formalism as above is Rabi spin oscillations in an- 
other guise; it utilizes NSI phenomenology which is con- 
tained implicitly in many published expositions. However 
we do not find another treatment which explicitly devel- 
ops the phenomenology of e^r to expressions which are 
as convenient as Eqs. (O, ([8]), and 

IV. IMPLICATIONS FOR CRUST TRAVERSAL 
A. MINOS baseline 

The recent MINOS results are based on NuMI neu- 
trino beam exposures of 1.7 x 10^" protons-on-target 



(90% CL) for survival. 

The apparent disparity between j/^ and data can 
be accounted for with the NSI matter effect. As seen 
from Eq. (O, the oscillation phase receives a matter- 
induced contribution which retards (advances) the to- 
tal phase for propagating (F^)- This shift in oscil- 
lation phase causes the observed survival probability to 
mimic vacuum oscillations with a Am^ of lower (higher) 
value. Consequently, the actual value of Am|2 accord- 
ing to our NSI formalism lies above (below) the best 
fit value of Ref. 0]. In the case of normal mass hi- 
erarchy, the NSI potential e^rVe needs to be negative 
(positive) for propagating (V^) and so e^r must be 
negatively-valued. On the other hand, for inverted hi- 
erarchy with Ato§2 < 0, C/iT would need to have a pos- 
itive value. For the MINOS 735 km baseline through 
the Earth's crust of density p ~ 2.72 g/cm'^, we use 
Ve ~ 1.1 X 10-" eV = (1/1900) km"^ 

For purposes of fitting, the slides of Ref. 01 were mag- 
nified to extract the data points and the error ranges 
presented. For instances of asymmetric errors, the larger 
of the two is chosen. A fit is performed over the MI- 
NOS data ratios which express the oscillation probabili- 
ties. The data is fitted to the phenomenological predic- 
tions of Eqs. 0, (HI), and ^ with Am^2, sin^ 2023, and 
(10) treated as free parameters. The role of these parame- 



ters is defined by the NSI phenomenology, which reduces 
to vacuum — >■ Vt mixing as {e^rl approaches zero. 



Figure 1(a) shows the MINOS neutrino data as ratios 
of the observed rate of events at 735 km to the Monte 
Carlo expectation in the absence of oscillations. The 
binned ratios are shown as a function of event recon- 
structed energy E^ for charged current (CC) events 
recorded in the MINOS far detector from FHC expo- 
sures totaling 7.2 x 10^° PoT. Superimposed (solid 
line) is the binned matter effect survival probability 
for the best-fit parameters as determined by both i/^ 



and data (described below), using Eq. 



We 



find: 



Am 



i2,sin2 20^^^ sR(e^,)|K|) 



(2.56 



+0.27 
-0.24 



10-3 eV^, 0.90 ± 0.05, -(0.12 ± 0.21) |V;|), where K is 
evaluated for terrestrial crust. Also shown for compari- 
son (dashed line) is the fit to z/^ — > vacuum oscilla- 
tions obtained by MINOS. The best fit parameters are 
(2.35 X 10-3 eV^ 1.0,0.0). 

Fig. 1(b) shows the corresponding distribution of ratios 
for antineutrino events observed over the null oscillation 
expectation, as a function of Ei, for j7^ charged current 
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FIG. 1. Histogram 1(a) (l(b)) shows the survival probability averaged per bin for i^fj, (77^ — s- 77^) at 735 km. The MINOS 

data points (Ref. 0]) show the ratios of observed CC-t'^ (CC-F^) event rates to rates predicted by MINOS in the absence of 
oscillations. In each Figure, the solid histogram shows the result of fitting the matter effect (77^) survival probability of Eq. 
(|8]) to the MINOS and oscillation survival ratios. The dashed histogram shows the vacuum oscillations (77^) survival 
probability for parameters of each of the MINOS fits (carried out separately for each sample). 



events of the MINOS RHC exposure of 1.7 x 10^" PoT 



The solid curve displays the best fit as shown in Fig. 1(a) 
differing only in the sign of e^r I I • Shown for compari- 
son (dashed hne) is the MINOS result obtained by fitting 
with 17^ — >• Vr vacuum oscillations; the best fit oscillation 
parameters are (3.36 x lO'^ eV^ 0.86, 0.0). 

Figures [2lja),[2{b) show confidence level (CL) contours 
for |e^T-| versus Am^ from fitting to the MINOS sam- 
ples with sin^ 29 marginalized (at each point, the non- 
displayed fit parameters are chosen so as to yield the 
lowest x^). In each plot, contours derived from the vari- 



ous fits are defined as the set of points on the |e^^| ver- 
sus Ato^ plane where Ax^, the difference between 
and the best fit is 2.3 (68% CL), 4.61 (90% CL), or 



9.21 (99% CL). Figure (Ha) shows |e^^| versus Am^ con- 
tours which are determined separately for the neutrino 
and antineutrino samples. Fitting to the data yields 
the ellipsoidal contours which surround the best-fit point 
(solid star) at (2.34 x 10~3gv2^ 0.94, -0.0041^1)- The 
contours exhibit a positive correlation between higher 
Am^ values and larger (more negative) |e^r|- Also 
shown are the corresponding contours obtained from fit- 
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FIG. 2. 

to the Vfi-CC and i/p-CC survival probability ratios as func- 
tions of observed at the MINOS far detector (Ref. 01). 
Plot (a) shows the contours from separate fits to (ellip- 
soidal regions) and 77,1 (shaded regions). Plot (b) shows the 
three-parameter combined fit to both and 77^. Here, the 
correlation properties exhibited by the separate fits of Plot (a) 
tend to merge for non-zero, negative values of e^r. In both 
plots the best fit points are denoted by stars; the up- (down-) 
pointing triangle in Plot (b) depicts the vacuum oscillation 
result of Ref. Hi for Uu (Va)- 



ting the z/^ data, which surround the best-fit point at 
(2.98 X 10-3eV^0.83,-0.23|ye|)■ In contrast to the cor- 
relation property of the contours, the contours 
correlate lower Am^ values with larger je^jrl- The 77^ 
contours are also much broader - less restrictive - than 
those from the i^^ fit, providing a striking visual reminder 
that event statistics are relatively poor compared to 
statistics in currently available long baseline expo- 
sures. Nevertheless one can see from Fig. [DJa) that the 
opposite-leaning correlations exhibited by the two inde- 
pendent sets of contours will tend to merge for non-zero, 
negative values of e^,- when both h'^ and 17^ are fitted 
simultaneously. 

Figure ^h) displays the Ax^ contours and the best 
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FIG. 3. The Ax for sin 2^23 with Am and e^r marginal- 
ized. The fit is to the MINOS data ratios versus Ei, which 
expresses i/^ and 77^ survival probabilities at 735 km (Ref. 2]). 
The solid (dashed) horizontal line shows the contour bound- 
aries at the 68% (90%) confidence level. 

termediate' Am^j value with the NSI matter coupling 
negatively- valued and of 12% of the strength of the con- 
ventional MSW matter effect. Figure [3] shows the Ax^ 
obtained from fitting to sin^ 2023 with the parameters of 
Figure HJb) marginalized. Non- maximal mixing, namely 
sin^ 2023 — 0.90, is preferred; sin^ 26*23 — 1-0 is excluded 
at greater than 90% CL. 

The per degrees of freedom for the three-parameter 
fit to both samples is 44.1/50 = 0.88. The fitted data 
comprises bin distributions of 1986 events together 
with only 97 of 17^ events. Obviously, the MINOS 
sample - at present - exerts dominant statistical power 
within the fit. Our three-parameter fit does not consti- 
tute supporting evidence for a claim concerning the 
NSI. On the contrary, it is readily discerned in Fig. ^h) 
that solutions with e^r = lie within the one-sigma con- 
tour, hence the NSI fit provides a description which is 
no better (or worse) that that afforded by conventional 
vacuum oscillations. Rather, the NSI fit to MINOS data 
serves to illustrate how the phenomenology may play out 
upon inclusion of larger 17^ oscillation samples. 

The dashed histograms of Figs. 1(a) and 1(b) show the 
descriptions afforded by the separate vacuum oscillation 
solutions obtained by the MINOS collaboration. Our NSI 
scenario provides a characterization of both data sets us- 
ing only three oscillation parameters whose relationship 
to each set is specified by the phenomenology proposed. 
This is in contrast to the vacuum oscillation descriptions 
using two different pairs of parameters, namely Am^, 
sin^ 26 and Am^, sin^ 29. 



B. K2K and T2K baselines 

In view of the e^r matter effect scenario presented here, 
it is sensible that the K2K experiment reported a value 



of Am^ higher than that obtained with MINOS 



fit location on the je^,-! versus Am plane from the joint That is, 
fit to the and 17^ data. The fit converges to an 'in- sin^ 26 = 



the K2K best fit value is 2.8 x 10"^ eV^ with 
1.0; at 90% CL, the allowed Am? range is be- 
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tween 1.9 and 3.8 x lO'^ eV^ with sin^ 26* > 0.6 
This is because the K2K basehne of 250 km involves less 
matter traversal than does neutrino propagation in MI- 
NOS, and also the neutrino beam spectrum is confined 
to lower energies, 0.5 < < 3 GeV. The result is that 
the vacuum part of the oscillation phase will dominate 
the perturbing term from the matter effect. In the fu- 
ture, when the T2K experiment with 295 km baseline 
obtains a sizable exposure, it will be of interest to in- 
clude those — >■ z/^ results into fitting with the NSI e^,- 
scenario. 

V. IMPLICATIONS FOR MANTLE 
TRAVERSAL 

A. Phenomenology 

The matter effect survival probabilities proposed for 
i^fi — >■ i^fi and for 77^ — >■ 77^ have implications for longer 
neutrino propagation baselines, such as those of atmo- 
spheric neutrinos which may travel thousands of kilome- 
ters through the terrestrial mantle and core. Given the 
elevated densities of the matter fields traversed, one in- 
fers from Eq. (O that apparent differences in the values 
of Am^ for j/^ versus oscillations may be more pro- 
nounced than for neutrino baselines confined to the ter- 
restrial crust. The effects described by Eqs. ([7]-|9l) can 
be elicited by considering separately the behaviors of the 
oscillation phase (p and of the probability factor T of Eq. 
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FIG. 4. The value of AM^g defined by Eq. ^ for propaga- 
tion in the mantle, as a function of neutrino energy. The NSI 
scenario suggests conventional analyses based upon vacuum 
oscillations may infer efi^ective mass-squared differences that 
differ for Vfj, and oscillations and exhibit dependence upon 
E^. The value Am32 = 2.6 x 10"^ eV^ is assumed for the 
plotted curves. 

To examine the consequences of the matter effect for 
Am^ values extracted by analyses which assume conven- 
tional oscillations, we construct, as a toy model, the ap- 
parent or effective AM^^. We equate the vacuum phase 
expression AAI^gL/4Ei, to the matter effect oscillation 
phase |7V|L of Eq. 0. Then, 



AM^ff = Amj^^l + sin 26* \Ve\ {AE^/Aral^) + 

+ {e^rVef ■ {AE,/Aml^)f (11) 

For ~ —0.12, the second term under the square root 
in Eq. (jlip is negative for and positive for . Equa- 
tion pT|) implies separate curves for the apparent AM"^^ 
as a function of i?,^, for and v^^ oscillations. The effect 
is illustrated in Fig. 4, which shows the AM'^^ curves 
for neutrino propagation through a terrestrial mantle of 
density 4.7 g/cm^. A treatment of atmospheric event 
samples for which no event-by-event separation of an- 
tineutrinos and neutrinos is made, will likely converge on 
an 'averaged' description with AAf^^g ~ Am^2- On the 
other hand, in investigations in which samples and 
i^f^ samples are separately analyzed, two or more AM^g 
values may be inferred, depending upon the E^ regimes 
which characterize the various event samples (e.g. sub- 
GeV, multi-GeV contained, multi-GeV partially con- 
tained, upward-stopping muons, upward through-going 
muons, etc). 

For neutrino propagation in matter there arises an ad- 
ditional complication due to the probability factor J- of 
Eq. For 623 mixing which differs from maximal, 

the J- factor introduces damping of the oscillation sin^ 0; 
the damping increases as cos^ 2^23 is moved away from 
0. Due to the sign difference in Eq. the damping is 
small for 77^ propagation. However, for z/^ propagation 
at energies for which the negative term within the de- 
nominator of Eq. ^ competes with the other two terms 
which are positive, the alteration of — ?■ survival 
probability can be dramatic. Figure [S] illustrates the os- 
cillation cloaking capability of the J- factor. For prop- 
agation in the terrestrial mantle, oscillations of i/^ with 
E,y Ri 29 GeV will be strongly damped with deviations 
from maximal mixing. This effect was identified and dis- 
cussed previously in Ref. [l^ whose authors designated 
it as an "anti-resonance" . 

The pivotal middle term within the denominator 
of Eq. ^ that gives rise to different damping be- 
havior of T for i/^ and cases is proportional to 
(e^rVe) X (^"^32)^^- We note that for the inverted 
hierarchy scenario which also fits the MINOS data, 
namely (Am§2,e,.T) = (-2.56 x 10-3eV^ -H0.12), the 
anti-resonance remains to affect the propagation, and 
so the factor e^rVe does not provide discrimination of the 
neutrino mass hierarchy. 

B. Atmospheric neutrino experiments 

The differences between and Vfj, propagation in ter- 
restrial mantle indicated by Figs. S] and [5j may mislead 
atmospheric neutrino analyses which use events having 
good resolution for L/Ei, but fit combined samples of 
(v^ +T'fj^) with vacuum i^^ v^- oscillation phenomenol- 
ogy. It is intriguing in this regard that the negative-log- 
likelihood contour for Avr? versus sin^ 1Q reported by 
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FIG. 5. Nine contour lines showing the magnitude of the 
probability factor F for i^^^ propagation in the mantle, Eq. 
(H]). The contours correspond to = 0.9, 0.8, 0.1. For 
atmospheric !/^'s with ~ 29 GeV and non-maximal 023 
mixing, the disappearance oscillation factor sin^ of Eq. ([9]) 
is strongly damped. 

the Soudan 2 collaboration does not exhibit a smoothly 
falling surface which converges to a single minimum. 
Rather, as shown in Fig. 9 of the experiment's final pub- 
lication 
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the contour exhibits two local minima in 
the regime near maximal mixing; one minimum lies above 
the current global-average 

Am^2 at 5.2 X 10"^ eV^, while 
the other falls below, at 1.7 x 10"'^ eV^. Similar behavior 
of a likelihood contour at maximal mixing was reported 
in the first MINOS analysis of 107 [v^ -\- 17, , ) contained- 
vertex atmospheric events (see Fig. 14 of [17l|). 

Clearly of interest are atmospheric event samples 
wherein 17^ and are distinguished, e.g. via magnetic 
tracking. The MINOS experiment has demonstrated 
such a capability with atmospheric neutrinos [Tsj . An 
analysis of a 24.6 kiloton-year exposure of the MINOS 
detector to atmospheric neutrinos is in progress 0. 

VI. DISCUSSION 

In view of the possible existence of an NSI matter ef- 
fect e^T- which may be large enough to distinguish 
from disappearance oscillations in current-generation 
long baseline experiments, certain near-term develop- 
ments are worthy of note: (z) As additional long baseline 

data from MINOS becomes available, more direct de- 
termination of limits or else measurement of the sign and 
magnitude of e^,- will become feasible. Currently, the MI- 
NOS collaboration is taking i/^ data with FHC running; 
plans are underway to enhance the 17^ sample with more 
RHC running, (ii) Analysis of atmospheric neutrinos 
in which F^'s are examined separately from f^'s using 

— > i^r vacuum oscillation phenomenology, may obtain 
Am^ which are larger than Am^ . If event samples repre- 
senting different energy regimes are separately analyzed, 
inferences of multiple A*m'^ values are possible, along 



the lines indicated by Fig. 4. (iii) An anti-resonance ef- 
fect is predicted for high energy i/^ oscillations in matter 
for the case of non-maximal 623 mixing. Its elucidation 
would be difficult for propagation baselines confined to 
the Earth's crust, however the effect might be accessible 
via study of i/^ propagation through the mantle. 

The efiT matter effect scenario of this work differs from 
the two other proposals which were motivated by the 
early MINOS results based upon FHC running. In 
CPT- violation hypotheses considered in Q, Am|2 and 
623 are to be distinguished from Am^2 ^^'^ ^23, whereas 
in our view such distinctions are a mirage arising from 
application of incomplete phenomenology. Our outlook 
is more nearly akin to the "apparent CPT violation" pro- 
posal of Ref. ^. However in the latter proposal, a light 
sterile neutrino conspires with a B — L interaction to 
generate a difference between I7^i and i^^ disappearance. 
A direct manifestation of mixing of active neutrinos to a 
sterile neutrino is that the ratio R of observed to expected 
neutral current rate at a far detector site is predicted to 
be less than 1.0 for either of or 77^ exposures In 
contrast, our NSI matter effect gives rise to probability 
oscillations among active neutrino fiavors, but with no 
net flavor loss, hence R measurements are expected to 
yield 1.0. 

Other proposals for new physics to distinguish from 
ly^ oscillations 23,, 2^ , and remarks on using NSI for this 
purpose [25], have appeared in recent preprints. Impli- 
cations for accelerator neutrino experiments arising from 
neutral current NSI, and a possible role for charged cur- 
rent NSI of the type v^-^N — ^ ^J■+X, are discussed in [26| . 
The phenomenological exposition presented here, as well 
as those to be found in all of the other above-mentioned 
works, have one theme in common: The data currently 
emerging from exposures at various baselines deserves 
careful scrutiny for evidence of particle interactions which 
lie outside the purview of the Standard Model. 
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